Superfluid density and competing orders in d-wave superconductors 
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We derive expressions for the superfluid density ps in the low-temperature limit T — + in d- 
wave superconductors, taking into account the presence of competing orders such as spin-density 
waves, idxy-psLinng, etc. Recent experimental data for the thermal conductivity and for elastic 
neutron scattering in La2-a;Sri:Cu04 suggest there are magnetic field induced anomalies that can 
be interpreted in terms of competing orders. We consider the implications of these results for the 
superfluid density and show in the case of competing spin-density wave order that the usual Volovik- 
like depletion of ps{H) is replaced by a slower dependence on applied magnetic field. We find 
that it is crucial to include the competing order parameter in the self-consistent equation for the 
impurity scattering rate. 

PACS numbers: 74.25.Nf, 74.72.Dn, 75.30.Fv 



I. INTRODUCTION 



The doping and temperature dependence of the superfluid density ps{x,T) and its correlation with the critical 
temperature Tc in high-temperature superconductors (HTSC) have been intensively studied since their discovery jj]. 
The last fewyears are also marked by new observations |M ^ 1^ particularly for the underdoped side of the phase 
I diagram This special attention payed to the superfluid density is not surprising, because all properties of 

'"^ , Ps{x,T) mentioned above show that HTSC differ significantly from previously known superconductors. An important 

■ example is the temperature dependence of the in-plane penetration depth A(r) which provided one of the first 
O \ strong evidence fo^' unconventional pairing symmetry and for the important role of gapless low-energy quasiparticle 

excitations. 

Another interesting dependence of ps is its variation with the external magnetic field H . This dependence contains 

■ information on the quasiparticles and on their interaction with the vortices present in the vortex phase of HTSC. 
^ ' Since these are extreme type-II superconductors, a huge vortex phase extends from the lower critical field, Hd ^ 

2^ . 10 — 100 Gauss to the upper critical field, Hc2 ~ 100 T. The leading term in the field dependence of Ps{H) ^ \fli is 
well understood Refs. IT^ using the local Doppler-shift approximation. This very useful approximation was 

introduced by Volovik |l5| . who predicted, that in contrast to conventional superconductors, in d-wave systems the 
density of states (DOS) is dominated by the contributions from the excited quasiparticle states rather than the bound 
\f~\ states associated with vortex cores. It was shown that in an applied magnetic field H the extended quasiparticles 
: DOS, 7V(w = 0, H) ^ ^/H rather than ^ H of the conventional case. The validity of the Doppler-shift approximation 
was discussed in Refs. 0, using the quasiclassical Eilenberger equations. It was shown Ref. ^3 that for the 
superfluid density this approximation works reasonably well at low temperatures. 

The characteristic behavior has been observed in the specific heat and thermal conductivity (see Ref. for 
a review of experiments). The dependence Ps{H) can also be directly extracted from mutual-inductance technique 
measurements fA Il9l l20| or, for example, from muon spin rotation (/iSR) measurements (see Ref. plj for a review). 
Q A big advantage of the former method is that it provides directly the desired dependence of ps{H). Moreover, the 
Q . measured inductance is directly related to the superfluid density, so that in contrast to the thermal conductivity 
^ ' discussed below there is no need to subtract a phonon contribution from the raw data |18| |. On the other hand, the 
• i-H , analysis of ^SR experiments involves additional assumptions about the internal magnetic field distribution in the 
' vortex state [2ll | which are contained in a model dependent function f{H). This important field dependent function 
H ] enters the relationship a{H) = f{H)ps{H) between the measured muon depolarization rate a{H) and ps{H). 
- - < Since 1997 there is ongoing discussion (see Ref. 122] for a historical overview) about the dependence of the thermal 



conductivity k{H) and its deviations from the expected ^ vH behavior. The crucial fact is that when the magnetic 
field is applied perpendicular to the Cu02 planes, the thermal conductivity shows a transition from a field-dependent 
regime k{H) ^ \Jli to a field- independent, plateau-like regime. The latest experimental results in underdoped 



'Electronic address: 'sharapo@mcmas ter.ca| 
t Electronic address: carbotte@mcmaster.cal 



2 



La2-:,Sr^Cu04 (LSCO) [H H at F = for < a; < 0.22 and [H HI at < F < 16T for 0.06 < a; < 0.22 are 
interpreted in Ref. in terms of a competing spin-density-wave order in the underdoped regime. 

A theoretical background for this interpretation was proposed in Ref. , where it was shown that in the presence 
of a spin-density wave gap m (2^ for = at T — > 



T " ~3 vpVA ^l + r 



(1) 



Here vp is the Fermi velocity, wa is the gap velocity and Tq is the impurity scattering rate which in Ref. is 
assumed to be independent of the gap m (we use units with h = c = 1, unless stated explicitly otherwise, and from 
Sec. Inl also set fcs = 1). 

As one can see from Eq. Q for m = and vp ^ va the minimal value of the thermal conductivity is Krnin/T = 
2kg /3. The opening of the gap m leads to a suppression of the thermal conductivity and allows for values of k/T which 
are less than Kmin/T, as is indeed observed in LSCO for x = 0.06 The presence of a nonzero field-dependent gap 
m{H) also allows one to explain the behavior of K{H,m{H)) observed in underdoped LSCO 25. 2(|. 

There are two additional experiments that support the idea of the presence of a competing order in LSCO. The first 
is an angle-resolved photoemission study Isoll that indicates the existence of a finite gap over the entire Brillouin zone, 
including the dx2_y2 nodal line in LSCO for x < 0.03. And an even more exciting observation was made for a LSCO 
sample with x = 0.144 and Tf, = 37K by elastic neutron scattering JSlI which showed that a static incommensurate 
spin-density- wave order develops above a critical field Hq « 3T. This picture is supported by the latest measurements 
[33 | made on a sample with the same doping which showed k{H)/T increasing for < H < 0.5T and decreasing at 
higher fields. However, a more heavily doped LSCO sample with x = 0.15 does not show any decrease of k{H)/T up 
to 17T. All these observations fuel interest in the investigation of competing orders. 

The purpose of the present work is to investigate the influence of competing orders on the superfluid density ps and 
its field dependence. We demonstrate that it is crucial to include the effect of the opening of the spin-density-wave 
gap m on the value of the impurity scattering rate T{m). When this effect is taken into account the the resulting 
dependence ps(H,m{H),r{m{H))) resembles the experimental results on very thin films of LSCO .20]. 

Although the latest experiments on LSCO [sil Is^ are mainly interpreted [23 in terms of a competing spin density 
wave (SDW) order, the earlier experiments on BSCO (see Refs. in and Ref. [s^l) were also interpreted using 
a complex d^y component generated in a d-wave superconductor in the magnetic field |3^ l35j . The removal of the 
d-wave node in optimally doped YBCO in a magnetic field was observed in the in-plane tunneling conductance |36j| . 
The authors of Ref. interpreted their observation in terms of competing idxy or is order. Here we also consider 
the influence of competing superconducting idxy and is orders on the superfluid density ps- This problem was 
considered previously by Modre et al. [3^ who calculated the London penetration depth for mixed superconducting 
order in zero magnetic field. We point out that the superfluid density ps may provide a way to distinguish SDW and 
superconducting orders. 

The paper is organized as follows. In Sec. ^ we introduce the 4x4 Dirac formalism convenient for the description of 
competing order with the underlying d-wave superconductivity. The general representation for the superfluid density 
is written in Sec. lIIII and in Sec. lIII Al the difference between ps{T — 0) for competing (with d-wave superconductivity) 
superconducting and SDW orders is considered. The dependence of the impurity scattering rate F on the values of 
competing gaps both in the Born and unitary limits is discussed in Sec. lIVI In Sec. we derive analytical expressions 
for Ps{T) in the presence of competing orders at 7J = in the low-temperature limit. The main results of the paper 
for the field dependence Ps{H) at T = are presented in Sec.^Jln Conclusions, Sec. I VIII we give a summary of the 
results obtained. 

II. DIRAC FORMALISM FOR DESCRIPTION OF COMPETING ORDERS IN d-WAVE 

SUPERCONDUCTORS 

We begin with the action for a d-wave superconductor written in imaginary time - momentum representation 

S = - j dT j d\<i^\T, k) [idr + T3^(k) - riA(k)] 7A(t, k), (2) 

where 

vl't(r,k)-(c|(r,k),ci(r,-k)) (3) 

is the Nambu spinor and cJ.(t, k) and Co-(t, k) with a =|, [ are, respectively, creation and annihilation operators. Most 
of the time we will rely on the nodal approximation, so that the precise form of the dispersion law of the quasiparticles 
^(k) and the d-wave superconducting gap A(k) is not essential. 
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It is impossible to consider other than is competing order while remaining within a 2 x 2 formalism. Depending 
on the physical assumptions made about the nature of the competing order , there are different possibilities for 
constructing a four-component field from Nambu spinors and switching to a 4 x 4 formalism (see e.g. jlH E3, ^3 ) ■ 
Since we are mostly interested in competing spin density wave which forms on top the superconducting state, we 
choose our spinors as was done in Rcfs. _28> .4L| 



^'l(t,k) = (4(t,k) cx(<,-k) c|(i,k-Q,) cx(i,-k + Q,) ) , (4) 

where = 2Ki is the wave vector that connects the nodes within the diagonal pair i = 1,2. Further since we 
are interested in the low-temperature (T <C Tc) properties of the system, we consider only the vicinity of the nodes 
k = Ki q with |q| < \Ki\ as shown in Fig.Hl Using that ^(k) = -^(k - Qj), and A(k) = -A(k - Q^) for k w K^, 
and then linearizing the spectrum as — vplx + 0{q^) and A(k) = VAly + 0{q^), one arrives at the low-energy 
action 

S = - j dr j dM!\{T, k) + Mivpqx + M2f Affy] *i(r, k) + (1 ^ 2, x < > y), (5) 

where Mi = as^r^ and M2 = — 0-3(8) ti, respectively. It is useful to reformulate the model (0) in the form of QED2+1, 
because this allows us to rely on the algebraic properties of 4 x 4 reducible representation of 7-matrices which satisfy 
Clifford (Dirac) algebra 

{7a<>7i^} = 2/45^^., = diag(l,-l,-l), = 0,1,2. (6) 

Another advantage of the QED2+1 formulation is that it also allows us to classify different competing orders in terms 
of different types of Dirac masses |2^ |3^ [s^ 0, . We introduce the Dirac conjugated spinor = ^-70, where 
7o is the 4x4 matrix that anticommutes with Mi and M2 and such that 7q =14. This choice is not unique, but we 
will follow the same conventions as in Refs. |23,|43 choosing 70 = cri (g) tq. Accordingly we define 71^2 via Mi = 7071 
and M2 — 7o72, so that 71 — ~ia2 T3 and 72 = i(J2 'Si ti satisfy Eq. J^. Finally we arrive at the action 

^ = - j dr J dk^i{T,k)[-fodr + jiVFQx + y2VAqy]'^i{T,k) + {1 — >2,x< — > y). (7) 

One may consider quasiparticle gaps rrii of different nature encoding it in the matrix structure Oi — 
(/4, 175, 73, 7375). Here the matrices 73 and 75, anticommuting with the matrices j^, and are 

73 = io-2 (E) T2, 75 = CT3 h- (8) 

Then, different gaps correspond to different types of Dirac masses added to the action In particular, the mass 
mi, with Oi = Ia, describes the (incommensurate) cos spin-density- wave (SDW), and the mass TO2 with O2 = «75, 
describes sin SDW. The masses my, and 7714 with O3 = 73 and O4 = 7375, correspond to the idxy-v^vc'vag and the 
js-pairing, respectively |28l l4l| . In the present paper we concentrate mainly on the SDW gap, mi = m with the 
corresponding bare Matsubara Green's function 



Go(z^„, k) = (9) 



Go(za;„, k) = . .2 ^ ,,2 .2 ^ ,,2 .2 ^ .2 ' 



iu^nlQ - Vpkljl ~ WAfc272 + ml 
UJn + v'^F^l + v\k2 + m 
and on idxy gap, TO3 = A^;^^ with the corresponding Green's function 

iUn' -IO - ypklll - f Afc272 - A rf^^73 

Finally we should define the electric current operator in 4 x 4 formalism. In Nambu formalism it reads 

j{T,ci = 0) = eJ dkvF{k)i;Hr,k)i2^{T,k), v^^(k)^^, dk^-0^, (11) 
and ^(k) — — ^(k — Qi), we arrive at the expression 



j(T,q = 0) = e / dkvF(k)*^(r,k)cr3 «) /2^'(r,k) 

JHBZ 

= e dkvF(k)4'(r,k)7o75^'(T-,k), 
Jhbz 



(12) 



'HBZ 

where the integration is over the halved Brillouin zone (HBZ), i.e., over the domain with ky > in Fig. ^ This 
implies that after the nodal approximation is used, one should include in the integration only two neighboring nodes 
with i = 1, 2 as reflected in Eq. {Tj), because the opposite nodes are already included in the 4x4 formalism. 
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III. GENERAL REPRESENTATION FOR THE SUPERFLUID DENSITY 

The superfluid stiffness (or density divided by the carrier mass m*) is given by |42| 

A:^"(T,i?)^^=r.,-A:f(r,iJ), (13) 

where is the diamagnetic (or stress) tensor and Ag is related to the London penetration depth Al by the standard 
expression Ag — c^/(47re^A|). In Eq. I|13|l A„ is the normal fluid density divided by the carrier mass, calculated 
within the "bubble approximation" with dressed fermion propagators (i.e., with self-energy S due to the scattering 
on impurities included) but neglecting vertex and Fermi liquid corrections, 

^^-T E / 7:7^^^F.(k)«F,(k)tr[G(zu;„,k)7o75G(»c.„,k)7o75]. (W) 

„JHBZ l^""] 



n— — oo 



As shown in |42l | the vertex corrections can be neglected if the impurity scattering potential is isotropic in k-space. 
Likewise the Fermi liquid corrections can be taken into account along the lines of Ref. ^3 ■ When S = it is more 
convenient to begin with the sum over Matsubara frequencies in Eq. (|14|) as done in Sec. IIII Al below. In a more 
generic case S ^ 0, the momentum integration has to be done first. Using the spectral representation for the Green's 
function G{iujn, k) with the spectral function given by the discontinuity of the fermion Green's function 

A{lu, k) = --^ [GR{iu + lO, k) -Ga{^- iO, k)] , (15) 
one can easily sum over Matsubara frequencies in Eq. (|14|l and represent A„ in the form [l^ 

A!.= / Tdc^tanh " 



HBZ (2^)^7-00 2T Am (16) 

X tr[GA(t^,k)7o75GA(w,k)7o75 - Gfl(w, k)7o75Gfl,(w, k)7o75] 

which is more convenient for the integration over k. In Eqs. IjlSI) and l|lt)|) G/j_^(cL',k) = G(ia;„ — > a; ± iO, k) are 
retarded and advanced Green's functions. 



A. Properties of As for different competing orders 

To simplify the formal consideration, in this section we will assume that there is perfect nesting, ^(k) — — ^(k— Q^) 
and A(k) — — A(k — Qi) for all values of the momentum k and not only in the vicinity of nodes. Then for the 
diamagnetic tensor we obtain 

d'^k i9^e(k) 



^ — ^tr[a3®r3G(iw„,k)7n] 



dkidkj 



HBz (2^)2 dhdkj E{k) 2T 



where depending on the kind of competing order i?(k) = . /^^(k) + A^(k) + m? + {AH }. Here and in what follows 

y 

the notation + {A^ } implies that either the competing SDW order with the gap ra or the competing idxy order 
with the gap /^a^^y is considered. 

For the second term of Eq. evaluating the trace in Eq. (|14|l we obtain 

One may notice that there is an important difference between SDW and id^y cases in the sign before m? and A^ in 
the numerator of the Matsubara sum. Evaluating the sum firstly for A| we arrive at 

^^1 7C^^i-^(k)^^F,(k)l i?(k)^w'e2(k) + A2(k)+A2 . (19) 

JHBZ l^^^j cosh -^fj^ * 
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Note that the same result holds for is order. Making an integration by parts one may check that the superfluid density 
Ps remains finite when A(k) = 0, because the second gap is also superconducting and — only when both gaps 
are zero, A(k) = ^d^y — 0. 
For the competing SDW order 

-^VF^{k)vFJ{k) 

It is easy to check that when A(k) = 0, while m remains finite, the superfluid density ps becomes zero. In contrast 
to the superfluid density, the thermal conductivity is blind with respect to quantum numbers distinguishing the gaps 
m and A^^^ (^s), so that Eq. ^ is valid for all competing SDW, is and idxy orders. 

Using the nodal approximation |42| one can estimate that the depletion of the condensate caused by developing 
SDW order at T = 

A„ = ^|m|«7V(0).|.M, (21) 

TTVA A 

where iV(0) is the density of states (DOS) per spin in the normal state and A is the amplitude of the d-wave gap. 



2m^ 



tanh 



£;(k) , c'(k) + A2(k) 1 



2T 



cosh^ J 



(20) 



IV. INFLUENCE OF THE SECONDARY GAP OPENING ON THE IMPURITY SCATTERING RATE 



In this section we consider the influence of nonmagnetic impurities on the residual scattering rate in the presence 
of a competing order. We begin with the Hamiltonian written in Nambu formalism 

H„„p- JdkJ dk'Vk,k'*^(T,k)T3«'(r,k'), (22) 

describing the interaction V{r) — J2i ^'5(r — rO with the positions of a random distribution of impurities. Accord- 
ingly in the 4x4 formalism 



H„ap^ f dk [ rfk'yk,k'«'^(T,k)/2<^T3^'(r,k') 
JHBZ JHBZ 

= / dk f dk'yk,k'*(r,k)(-757i)*(r,k'). 
Jhbz Jhbz 



(23) 



/HBZ J HBZ 

Therefore the set of equations for the impurity self-energy reads 

Go\tu, k) = lujo - ^(k)7i - A(k)72 - mh - {Ad.„73}, (24a) 
G-\iCj,k) =Go\iCj,k)-Y.{iiJj), (24b) 
E(ia)) = r™p(-757i)[c + 5(ia))757i]"\ (24c) 

.(-) = ^/(0G(...,k), (24d) 

where c — 1/ {'kN{0)V) is the parameter which controls the strength of impurity scattering and Timp — nimp / {tt N (0)) 
with riimp being the concentration of impurities. Since here we used an unexpanded dispersion f (k) and d-wave gap 
A(k), the integration in Eq. Ij24d|l is done over the whole Brillouin zone to count correctly the contribution from all 
nodes. 

Expanding the self-energy E and g{iu}) in the 7-matrices 

E(icD) = E"(icD)7o + ^^{iCj)-f2 + ^\iCj)h + {I]3(iw)73}, 

(25) 

g{iLb) = 5"(itI')7o -l- g^(j'I')72 + 9 {iu})h + {9^ {i^jha} , 
one can obtain the system of self-consistent equations for 

iu) = iuj — TP{iLj), 
A(k,iw) = A{k) + J:^{iCj), 

m(ia)) = m + I]"^(ia)), (26) 

|A,^„(k,zcD) = Ad^^(k) + E3(,^)|. 
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Assuming particle-hole symmetry for simplicity, the renormalization of ^ is zero. In particular, when the competing 
orders are absent, the system of equations reduces to T-matrix equations for d-wave superconductor studied, for 
example, in Ref. j43|- Since in this case the averaging over the Fermi surface gives g^{iLo) = 0, the only relevant 
equation left is for a) or E". 

Here our goal is to take into account the influence of a competing order on Co. When a competing order develops, it 
affects both the equation mentioned above for cj [YP) and the new equation for Yj^ {iCj) {I]'^(ia))}. However, because we 
do have an explicit gap equation for m {lS.d^y} (see e.g. Ref. ^3); what follows we will assume that the dependence 
m,{ui{uj = 0)) is given phenomenologically and do not consider an equation for Y.^ (iui) {E^(ia))}. This assumption 
means that we do not distinguish the values m(a)(cL; ~ 0)) and rri, so that in what follows we denote the competing 
gaps as TO, Ad^^ and A.,. These gaps already include the effects of impurities and magnetic field and correspond to 
their phenomenological values extracted from experiment. Although in what follows we do not consider the equation 
for it is useful to stress the analogy between the the effect of nonmagnetic impurities on SDW order, and 

magnetic impurities in conventional s-wave superconductors that lead to the finite density of states inside the 
gap. Physically this means that the scattering off random impurities prefers to make the system homogeneous by 
washing out the nonuniform SDW structure |4q . 

Two other important assumptions that we make are the following: 

(i) We are interested in the value of the zero energy impurity scattering rate, F = — ImS5j(a)(u; — > 0)) and do not 
consider the effects related to the energy dependence of the self-energy E'^((D(w)) studied in Ref. [43. This assumption 
is justified by the fact that LSCO compound is intrinsically more dirty system than YBCO considered in Ref. 47] . 

(ii) Although we assume that the values of the competing gaps are field dependent, e.g. to = m{H) we do not include 
the inffuence of the Doppler shift on F ^3, Oil ■ 



A. T-matrix equation for competing SDW order 



We begin with the equation for TP for competing SDW order 



- imp 



1 



1 



1 



gO{Cu) + gi{Cj) 
I 



g'^{Cj)+gi{Co) c + gO{Cu)+gi{u) 



(27) 



which is obtained from Eq. (|24c|) with g{iui — ^ uj) given by Eq. (j24d|l using the decompositions H2t)|l . 



Since we are interested in the value S°(ti)) with a) 



E (J)) at w = 0, we need only the functions g (uj) and 



g {Co) calculated for uj = iV. Using the nodal approximation we obtain 



5°(*r) = 



1 



1 



7riV(0) ttvfVa 

1 I 
TrN{0) ttvfva 



In- 



Pi 



-iT, 



In- 



(28a) 
(28b) 



where the ultraviolet cutoff po is introduced. Note that for po ^ 4A and (vfva)^^ ^ TrN{0)/ A, where A is the am- 
plitude of the d-wave gap A(k) = A/2(cos k^a — cos fc„a), the function g{Lu) calculated using the nodal approximation 
agrees with the expressions given in Refs. ^3, 47, 48]. In Eq. H27p the Born limit corresponds to ^ (c ^ 00), 
while the unitary limit corresponds to — > 00, i.e. c 0. 



1. Bom limit 



In the Born limit Eq. (|27|l reduces to 



-V(-). 



(29) 



Due to the fact that only 5° enters Eq. H29|) , the Born limit appears to be the same both for SDW and superconducting 
orders, because g^ depends either on m? or on A^ . 

Substituting Eq. (|28a|) in Eq. (|29|) and solving it with respect to F we obtain 



- Po exp 



F 



'imp 



(30) 
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For a given impurity concentration this solution is nonzero only if m < rricr, where rricr = 
Po exp[— c^7r^?;^tiA-^(0)/(2rimp)]. Because the Born limit is considered, the value of rricr is exponentially small. 
Nevertheless in the unitary limit there is a possibility to have both finite T and rather large values of m. 



2. Unitary limit 

LSCO compound is intrinsically a dirtier s^tem than other cuprates, so that the unitary limit is more relevant. 
Moreover, thin LSCO films studied in Ref. [23| seems to have particularly large values of F ~ 6 -f- 50K which also 
indicates the relevance of the unitary limit. When c — > Eq. (|27|l reduces to 



^liQ) - F™p j^.^.^j^ [ [g^{Co)Y ■ ^^^^ 



Now because still enters Eq. the resulting transcendental equation for F is 



F' = T,^pT:'N{Q)vFVA 



In. ^0 



m\ (32) 



F^ + 77! 

which differs from that of Ref. by the last term m?. Due to its presence the dependence F(m) is rather strong 
(see Fig. 12) and as in the Born limit there is a critical value m™'* such that F(to) = for m > m™*. As we will 
argue below the observed deviations of the dependence Pb{H) (or As{H)) from ~ in high fields can be caused by 
the fact that the field dependence of m{H) affects the behavior of Ps{H) via the dependence T{m{H)). Finally we 
note that the case of competing is order is also described by Eq. H32|l . 

B. T-matrix equation for competing id^y order 

Similarly to Eq. H27|) for the competing id,j.y order we obtain 

d - [5"(t^)]^ + W{^)Y 

where g^{iT) is given by (|28a|l . while for dxy order g'^ — 0, analogously to the function g^ for dx2_y2 case. 



1. Born limit 



As was already mentioned in Sec. IIV A II in the Born limit there is no difference between the consideration of 
competing SDW and idxy order, so that one may simply replace the gap m by A^^^ in the corresponding equations. 
Moreover, this consideration is also valid for a competing is order. Competing is order with d-wave superconduc- 
tivity was considered in Ref. [s^l , where besides the T-matrix equations for the impurity scattering rate, the optical 
conductivity order was studied. 



2. Unitary limit 
Since g^ — 0, instead of Eq. we arrive at the equation ^ 



In 



(34) 



Its solution T{Ad.^^) is shown in Fig. [51 It demonstrates that a competing idxv does not significantly perturb the 
value of F with respect to the A^^ = case, when it reduces to the form 



F = -"^inipVEVA 



In^ 

F. 



(35) 
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SUPERFLUID DENSITY FOR H ^0 



In Sec. IIIlXl we discussed in a more formal way how the different competing orders affect the superffuid density. 
Here instead we concentrate on the simple analytical expressions that demonstrate the dependence of A„ on the 
impurity scattering rate T and the competing gaps m and A^j^ . In the nodal approximation the representation Hlt)|) 
acquires a form convenient for analytical calculations 



where 



with 



and 



A --^^J 

ZTTVA 



J = — doj tanh — -I (uj) 
Jo 2T 



1 fP" 
2vri Jo 



(36) 



(37) 



(38) 



(39) 



Here p = ^/ v'j^k^ + y-j^ is the dispersion law of the quasiparticles in the nodal approximation. J in Eq. H36|l is twice 
bigger than in Rcf. because a 4 x 4 formalism is used. Now we consider the cases of competing SDW, id^y and 
is orders. 



A. Competing SDW order 

Substituting the Green's function with the self-energy ImUfj Ai'^ = 0) = in Eq. (|39(l one obtains 



Ir,a{(^,p) 



Then integration over the energy p we arrive at 



4[{lu±iT)^+P^ 

[(W±ir)2-p2, 



arctan 



pI + t' 



2rw 



arctan ■ 



pgr 



2rw 



Finally integrating over lo for F 3> T (also for pq ^ m, F) we obtain 



A„ 



VF 
TT^VA 



Fin- 



pI 



F2 



m arctan 



2mT 



F2 



7r|m|0(m^ 



y2p 



3 F2 



(40) 



(41) 



(42) 



where 9 is the step function. The first term of Eq. (|42|l can be interpreted as the DOS contribution to the depletion 
of the condensate even at T = 0, because the DOS (per spin) in a dirty d-wave superconductor with a finite gap m 
reads ,2^] 



N„,{0) = 



-Fin- 



Pa 



Vr2 



(43) 



The second and third terms of Eq. H42() describe the depletion of the condensate because of the development of SDW 
order and for F — > they reduce to Eq. H21|l discussed above. Finally the last term of Eq. (|42|l ^ shows that the 
characteristic for d-wave superconductors i.e. a linear T dependence for A„ changes in the presence of impurities and 
becomes ~ T^. This behavior is indeed observed in thin films 0| over a wide range of the temperatures. Based on 
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Eq. H44|l below, in some films it was estimated that F can be as big as ~ 46. 6K '2^1 . When m 
to the known [ij, |52, |53 expression 



A„(r) 



Vf 



2r .1- /TJ 
TT r 3 r V r3 



0, Eq. (gSJ reduces 
(44) 



Again the first term of Eq. H44I) is proportional to the DOS of a dirty d-wave superconductor |43 |. 

It is also instructive to consider the limit F ^ when the T-dependence of A„ is expected to change from quadratic 
to linear. To extract this limit one should extract the singular part of /(w) in the limit F ^ , viz. 



i{LL>) = -sgnw[sgn(po - t^^) - sgn(m^ - w^)] + 2po[S{ll> - po) - S{uj + po) 
where we again took pq ^ m, F. Then 



J(ra) 



-4rin 2cosh- 



2^) 



Obviously for m = we recover the well-known expression 

An = -T. 



TT VA 



For \m\ ^ T we obtain 



A„ 



VF 
TTVA ' 



[|m| + 2rexp(-|m|/r)], 



where the first term coincides with Eq. (|21|l . For |77i| <C T, we arrive at the expression 



A, 



TT VA 



T 



1 



1 

\W2J^ 



(45) 
(46) 

(47) 

(48) 

(49) 



Here we refer for comparison to simple expressions for A„ when c?-wave superconducting states coexists with the 
orbital antiferromagnetic (d-density-wave) state [53 |. 



B. Competing idxy and is orders 

The difference between competing SDW and superconducting orders can be traced back to the opposite sign before 
the m? term in the numerator of Eq. H4U|) and the corresponding sign before in the formula 



/k,a(w,_p) 



4[(a;±zF)^+ p^ + Ag^ 



(50) 



which is related to the different 7-matrix before A^^ in Eq. (|10|l . Note that for competing is order, Eq. H5l)|l turns out 
to be the same, so that these order differ by the influence of the superconducting gaps A^j^^ and A^ on the impurity 
scattering rate F (see Sec. IIV K% in the unitary limit. Integrating over p one obtains 



/(a;)= — 

TT 



arctan ■ 



- UJ 



A2 



2Fw 



arctan ■ 



2Fw 



FA. 



(51) 



^pI + Al„)^ + r2 (^ + ^pI + AIJ^ + F2 
Finally integrating over lu for F » T (also for pQ Ad^,,, F) we obtain 



A„ 



VF 
TT-^VA 



Fin 



Pi 



r 



Al 



F2 



F2 



AL (r^ + AL)' 



(52) 
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Similarly to Eq. in the limit F ^ from Eq. (|FT|l we obtain 

i{uj) = - sgnw[sgn(p2 _ uj^) _ sgn(A^ - u^)] 



(53) 

+ 2po[S{uj - po) - 5{uj + po)] - 2Ad^^ [S{lj ~ A^^ J - S{uj + A^^ J], 
where we again took po ^ A^^^^, F. Then 

J{A^J = 2A,^^ tanh ^ - 4Tln (^2 cosh , (54) 

which differs from Eq. (|46|l by the term with tanh. For A^^,, = we again recover Eq. H47|l . For A^^^^ ^ T the 
dependence A„(r) becomes thermally activated due to the secondary gap 

2 vj^ 

An = -— Ad,„ exp(-Arf^„/r), (55) 
while for A^-^^ ^ T the leading term of A„(r) coincides with Eq. (|47|l 

1 



21n2wF 
A„ -T 



1 - 



(56) 



8 In 2 T2 

Comparing the last equation with Eq. H49|) one may notice that they differ by the sign before the term. 

VI. SUPERFLUID DENSITY IN THE VORTEX STATE 

The presence of circulating supercurrent around vortices in the mixed state can be taken into account in the 
semiclassical approach by introducing the Doppler shift in quasiparticle energies, uj ^ ui — Vs(r)k Here Vs(r) is 
the superfluid velocity at a position r which depends on the form of the vortex distribution and k is the quasiparticle 
momentum which can be approximated by its value at the node [5ll | . This distribution is described by the function 

'^(') = jj d'rS{e -Mr)k), (57) 

where the integration is over the unit vortex cell with the area A = ttR^. Several choices for Vie) were discussed in 
Ref. [SJ ■ Among them are the distribution for the vortex liquid [HJl 



which is the most convenient for analytic calculations, and |51ll5 



1 / e= 



for the completely disordered vortex state. The characteristic energy scale Eh in Eqs. (|58|l and (|59|l is associated 
with the Doppler shift energy in the vortex state 

Eu{H)=a— = a—^ — ^ (60) 

where a is a geometrical factor of order unity and H is the magnetic field applied perpendicularly to the ab plane. In 
the second equality we used the convention of Ref. [sJl (see also Ref. [5i| and Ref. [33] therein) that for a = 1 there is 
one flux quantum $o = hc/2e per unit cell of the vortex lattice approximated by a circle of radius R = {^q/ttHY/'^ . 
The final results depend somewhat on the choice of the distribution function and on the value of a, however, the 
qualitative results are not sensitive to this choice. In what follows we take the value vp ~ 2.5 x lO^cm/s '2^ which 
corresponds to Eh[K] = 38K • T-^/^yi/fr] as used in Ref. 
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Now the Doppler shift effect can be incorporated in the Green's function formaUsm that already includes the 
scattering on impurities by averaging over the distribution V{€) [5llf57| 



(61) 

X tr[GA(t^ - e,k)7o75GA('^ - e,k)7o75 - Gr{uj - e, k)7o75G'fl,(w - e,k)7o75]. 

We note that in 2 x 2 Nanibu formalism the replacement to ^ lo — e for taking into account the Doppler shift in 
the argument of the Green's function is exact. In 4 x 4 formalism the corresponding Green's function describes two 
opposite nodes with the reversed sign of the Doppler shift. Nevertheless one can check that in the latter case the 
simple prescription uj ^ uj — e \s approximately valid if we neglect the terms ~ em. In the nodal approximation 
similarly to Eq. (|36() one may write 

/>oo 

A„(i/) = --^/ deV{t)J{m{H)-e), (62) 

where the function J(e; m{H)) is given below. However, the field dependence of A„(iJ) is not yet completely specified, 
because one should also provide a field dependence for the gap m{H) {A^^ (H)}. We will come to this question later 
in Sec. Em 



A. Results for the field dependence of the superfluid density for m, A^^ , F = const 



For T = and pq ^ to, Ad^^ , F in SDW case one obtains 

^2 ^2 



J(m,T ^ 0;e) 



e arctan 
+Fln 



2eF 



2l^l 



+m arctan 



V[F2 + (TO-e)2][r2 + (TO + e)2] 
2mF 



(63) 



F2 



'K\m\9{m 



and in the id^y case 



J(Arf^„,T = 0;e) 



e arctan ■ 



2eF 



+Fln- 



V[F2 + (A,^„-e)2J[r2 + (A,^„+6) 



(64) 



Before doing the numerical calculation it is useful to consider analytically the limit F — > 0. One may notice that 
in this limit the only difference between the SDW and idxy cases is the last two terms to in Eq. (|63|l . Thus we 
consider firstly the idxy case and then discuss the role of these terms. We obtain from Eq. ()64|l that 



Substituting Eq. ^ in Eq. ^ and using the vortex liquid distribution H58|) we arrive at 



El 



(65) 



(66) 



This means that in spite of the presence of a gap for nodal quasiparticles, the contribution to A„ still has the behavior 
~ ^/H if Ad^y (H) < ^/H. The origin of the gapless behavior is the Doppler shift of the quasiparticle energy i?— kv^ (r) , 
which is position dependent. There are regions where the shift is larger than the minimal gap Ad^y in the spectrum, 
thus leading to the finite DOS and A„. This point about the DOS was emphasized in Refs. [sJIs^, where it is stressed 
that regardless of the power with which Ad^^ ~ opens up in the field, as long as a > 1/2, the leading term in 
DOS will always be ^/H at small fields. Thus in the clean system A„(7J) would remain ^ ^/H even in the presence 
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of a nonzero gap A^^j, oc H" with a > 1/2. As we already mentioned in the SDW case the last two terms of Eq. H63|) 
contribute in A„. However, if the SDW gap m{H) cc V^, the Vll behavior of A„(iJ) would persist, because the 
contribution of the last two terms of Eq. is vp / {T^Vi^)w?{H) / E'jj + m?{H). 

In Figs. 131 and 21 we plot the dependence As{H) for constant field independent gaps m, A^^^ for the clean case 
(r = 0) and for T = 16K, respectively. This dependence is obtained numerically from Eq. H62() using the distribution 
Eq. H59|l . We use vf/va = 30 and a rough estimate of the diamagnetic term r = 1500K from the Uemura plot. 
In the clean limit F = at zero temperature, the opening of the gap leaves the condensate unaltered if it is of 
superconducting character, but depletes it for the spin density wave case. The application of an external magnetic 
field H oriented perpendicular to the cooper oxide planes creates quasiparticles and decreases the condensate density 
as compared with its zero field value. Comparing with the pure dr^-2_y-2 case, the reduction in condensate for a given 
value of H is less for an additional id^y (is) gap and even less for the SDW case. However, because in this last case 
the H = value is already depleted as compared with the pure dx2_y2 case, the condensate density remains lower for 
all values of H considered. Looking at Fig. (31 one may develop the impression that the cases of competing id^y and 
SDW orders are easily distinguishable experimentally. However, the situation is more complicated because in practice 
one considers the dependence As{H) — As{H = 0) which can hardly distinguish different competing orders. Moreover, 
in all cases the introduction of impurity (residual) scattering reduces the condensate density and makes the effect of 
H on it more similar (Fig. . 

B. Ansatz for gap ni{H) 

In general the competing gap and its doping dependence have to be obtained by solving a self-consistent system of 
equations for the main d-wave and subdominant gaps (see e.g. Ref. 1^ H 1 . However . since the origin of the secondary 
gap is unknown, here we follow the phenomenological approach of Ref. where it was assumed that 

m{H, a;) = (1 - a;/0.16)i/26'(0.16 - x){mQ + bEn), (67) 

where x is the doping, mo and b are free parameters. 

It was demonstrated in Ref. [2^ that the experimental data in Refs. H^H^ can be qualitatively understood using 
the gap l|H7|) which is generated below a critical doping Xcr = 0.16 and increases with inagnetic field as VH [131 ■ As 
mentioned in the introduction, the present work is partly motivated by the experiments [8lll82|| made on an a; = 0.144 
LSCO sample, where the SDW gap develops above a critical field Hq « ST. Since here we do not consider the doping 
dependence of the gap but instead are more interested in the effects related to the critical field, we assume that 

m{H) = bEH{H - Ho). (68) 

Moreover, while in Ref.j22| a rather large value of 5 = 2.2 was used, in the present paper we consider the case of 
small values of 6 = 0.17, so that m{H) <C Eh- The dependence m(H) is shown in Fig. |S1 where we also plot the 
dependence T{H) obtained by solving Eq. for r(m) with po = 250K and n'^TimpN{0)vFVA = 1500K^. As we 
already mentioned, there is also a direct influence of the Doppler shift and Andreev scattering on F ■ As 

shown in Ref. in the unitary limit the change of F due to the Doppler shift is not significant. Nevertheless these 
effects and particularly the energy dependence of T{H,ui) j4^ will become important when the value of to < m™* 
and r{uj — 0) approaches zero. 

Finally we note the fact that the explanation of experimental data requires a doping dependent gap, which supports 
its SDW character. The generation of a idxy gap by a magnetic field can presumably occur at any doping. 

C. Results for the field dependence of the superfluid density for field dependent gaps and F 

Substituting into m{H) and T{H) shown in Fig. O for the SDW case and into Ad^^^iH) = m{H) and 
F = const = T{/S.d^^{H = 0)) for the idxy case and using Eq. Ht)2|l after numerical integration over e with V{e) 
given by Eq. (|59|l we obtain the results shown in Fig. |B1 We use vf/va = 30 and r as above. Above H = 3T the 
difference between Ps{H) without competing order and Ps{H) for the idxy order is hardly noticeable, because we 
used a small value of 6 = 0.17. Nevertheless, for the SDW order Pa{H) in high-fields deviates from a i/ff behavior 
quite significantly. This effect is caused by the decrease of F(i?) seen in Fig. O We stress that a similar behavior of 
As {H) — As {H = 0) with a crossover from y/H behavior in low fields to In H dependence in high fields is observed in 
one of the samples with the lowest Tc studied in Ref. |23|. 

In contrast to Ref. 28], where the decrease of the thermal conductivity k{H) is directly associated with the opening 
of the gap and, accordingly, requires rather large values of b, the results presented here for Ps{H) are related to the 
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indirect influence of the development of the gap on T{m{H)) which in turn leads to the deviations of the dependence 
Pn{H) from a simple ^ \/H law. 

VII. CONCLUSIONS 

We have addressed the problem of the superfluid density of a d-wave superconductor with competing order. Both the 
case of a spin density wave (SDW) and a second (minority) superconducting order with idxy symmetry are treated and 
compared. The nodal approximation is introduced to treat the main d-wave gap and so the formulation is restricted 
to low temperatures. The resulting action which corresponds to QED2+1 involves a reducible 4x4 representation of 
Dirac matrices with competing order equivalent to the different types of Dirac masses. Residual impurity scattering 
is accounted for within a T-matrix formalism which includes as special cases the Born and unitary limit. 

SDW and a second superconducting order are, in principle, very different. For example, in the limit when the 
main d-wave order is set zero, the superfluid density vanishes for the SDW case, but remains finite for the idxy 
(or is) superconducting case. Also it is found that the SDW order reduces the superfluid density at T = (see 
Eq. 1)2111 ') because it competes for fermi surface with the d-wave order, but for idxy order there is no such effect. At 
low temperatures T <^ m, A^^ the superfluid density acquires an exponentially activated form (cp. Eqs. (|48|) and 
(|55|l 'l. In the opposite limit T <C m, A^^^ <C A, where A is the amplitude of the main d-wave gap, the well-known 
linear in T law (|47|) is modifled (see Eqs. ()49|l and (|55|l ) and an additional 1/T dependence appears with a coefficient 
proportional to the square of the secondary gap value. The sign of the correction is opposite in the two cases. So far 
we described results only for the pure limit. Analytic results are also obtained for the case when impurity scattering 
is in the limit when the zero frequency value of the impurity self-energy F 3> T. The expressions properly reduce to 
the known results when the secondary gap is set to zero. In the both cases impurities modify the classic linear in T 
dependence to a dependence as is also the case in d-wave superconductor without competing order. The coefficient 
of term (see Eqs. (|42f) and (|52|l 'l however is modified by the presence of the secondary gap and this modification is 
different for SDW and idxy superconducting order. This is also the case for the change of the zero temperature limit 
of the superfluid density due to impurity scattering. 

We have also considered the influence of the opening of a secondary gap on the magnetic fleld dependence of the 
superfluid density for H oriented perpendicular to the Cu02 plane. We have found that the presence of competing 
orders causes deviations of the fleld dependence of the superfluid density Ps{H) from a simple law which is 
associated with the Doppler shift of the quasiparticle energy. A new and rather important conclusion is that this 
effect is caused not only by the developing competing order, but also by the influence of this order on the residual 
impurity scatte ring rate at zero frequency. On the experimental side the superfluid density is a directly measurable 
quantity 0)E)l23, so that the effects discussed here may be well within experimental access. Another advantage 
of the superfluid density is that it is not blind with regard to the different kinds of competing orders as is thermal 
transport. However, on a theoretical side, in contrast to the thermal transport, charge transport is renormalized by 
the vertex and Fermi liquid corrections |42l | which were not considered in the present paper. Accordingly the values of 
vf^ i'A) F and m extracted from the measurements of the superfluid density may be in disagreement with the values 
for the same parameters extracted, for example, from the measurements of thermal conductivity. 
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FIG. 1; (Color online) The schematic graph of the Fermi surface with the vectors Ki, i — 1,2, and q. The SDW ordering 
vectors are Qi — 2Ki. 
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FIG. 2: (Color online) The dependence of the impurity scattering rate F on the value of the gaps m and Ad^y of the developing 
SDW and idxy orders, respectively. For m > rrict the value F(m) = 0. 
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FIG. 3: (Color online) The dependence As{H))/t at T = for d-wave superconductor with additional SDW, id^y gaps 
^dxH = ra = 15K and also without any competing order in the clean, V = OK limit. 
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FIG. 4: (Color online) The dependence As{H))/t at T = for d-wave superconductor with additional SDW, idxy gaps 
^dxM = m = 15K and also without any competing order. The constant F = 16K is taken. 



17 




FIG. 5: (Color online) The model dependence of the gap m on the applied field H and the resulting dependence of the impurity 
scattering rate F on m. 
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FIG. 6: (Color online) The dependence (As(0) — Aa{H))/T at T = for d-wave superconductor with additional SDW, id,, 
gaps (see Eq. I|68|l and Fig. |3 and also without any competing order. 



